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^ , Abstract 

We construct XX- and Hubbard- like models based on unitary superalgebras gl{N\M) 
^ ■ generalising Shastry's and Maassarani's approach of the algebraic case. We introduce 

the R-matrix of the gl{N\M) XX model and that of the Hubbard model defined by 
coupling two independent XX models. In both cases, we show that the R-matrices 
satisfy the Yang-Baxter equation, we derive the corresponding local Hamiltonian in 
the transfer matrix formalism and we determine the symmetry of the Hamiltonian. 
Explicit examples are worked out. In the cases of the gl{l\2) and gl(2\2) Hubbard 
models, a perturbative calculation at two loops d la Klein and Seitz is performed. 
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1 Introduction 



The Hubbard model was introduced in order to study strongly correlated electrons [1,2] and 
has been used to describe the Mott metal-insulator transition [3,4], high T c superconductivity 
[5,6] and chemical properties of aromatic molecules [7]. Since then, it has been widely 
studied, essentially due to its connection with condensed matter physics. The literature on 
the Hubbard model being rather large, we do not aim at being exhaustive and rather refer 
to the books [8, 9] and references therein. Exact results have been mostly obtained in the 
case of the one- dimensional model, which enters the framework of our study. In particular, 
the ID model has been solved by means of the Bethe ansatz in the celebrated paper by 
Lieb and Wu [10]. However, the set of eigenfunctions considered there was incomplete, and 
a complete set of eigenstates was constructed in [11] using the 5*0(4) symmetry of the ID 
Hubbard Hamiltonian. 

Although the Hubbard model certainly exhibits fascinating features among integrable 
systems, the understanding of the model within the framework of the quantum inverse 
scattering method appeared only in the mid eighties. The R-matrix of the Hubbard model 
was first constructed by Shastry [12,13] and Olmedilla et al. [14], by coupling (decorated) R- 
matrices of two independent XX models, through a term depending on the coupling constant 
U of the Hubbard potential. The proof of the Yang-Baxter relation for the corresponding 
R-matrix was given by Shiroishi and Wadati [15]. The construction of the R-matrix was 
then generalised in the gl(N) case by Maassarani et al., first for the XX model [16] and then 
for the gl(N) Hubbard model [17,18]. Within the QISM framework, the eigenvalues of the 
transfer matrix of the Hubbard model were found using the algebraic Bethe ansatz together 
with certain analytic properties in [19-21]. 

One of the main motivations for the present study of the Hubbard model and its gener- 
alisations is the fact that it has recently appeared in the context of N = 4 super Yang-Mills 
theory in two distinct ways. Firstly, it was noticed in [22] that the Hubbard model at half- 
filling, when treated perturbatively in the coupling, reproduces the long-ranged integrable 
spin chain of [23] as an effective theory. It thus provides a localisation of the long-ranged spin 
chain model and gives a potential solution to the problem of describing interactions which 
are longer than the length of the spin chain. The Hamiltonian of this chain was conjectured 
in [23] to be an all-order description of the dilatation operator of N = 4 super Yang-Mills 
in the su(2) subsector. That is, the energies of the spin chain are conjectured to be the 
anomalous dimensions of the gauge theory operators in this subsector. In relation to this, 
an interesting approach to the Hubbard model is given in [24] that leads to the evaluation 
of energies for the antiferromagnetic state and allows one to control the order of the limits 
of large coupling and large length of the operators/large angular momentum. 

The Hubbard model has also arisen in a slightly different way in the context of N = 4 su- 
per Yang-Mills (SYM). Following reasoning developed in [25], the long range spin chain 
describing N = 4 super Yang-Mills theory can be described in terms of scattering of 
momentum-carrying excitations (at least in the limit of very long operators or chains). Under 
the assumption of integrability, this scattering is governed by a two particle scattering ma- 
trix which is essentially determined up to an overall phase factor by s-u(2|2) symmetry [26]. 
This phase factor was introduced in [27] where its importance for matching with data from 
the string theory regime was discussed. In a recent paper [28] it has been shown that the 
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S-matrix thus derived satisfies the Yang-Baxter relation (or a twisted version, see [29]) and 
in fact is proportional to the tensor product of two copies of Shastry's R-matrix [12, 13]. The 
undetermined dressing phase of the S-matrix can be constrained by appealing to crossing- 
symmetry [30]. A proposal for its complete form was given in [31] based on an earlier 
guess [32] and conjectures for the form of the string Bethe ansatz [33]. The non-triviality of 
the dressing phase leads to modifications of the proposal of [23] at four loops and beyond. 
Following the suggestion of [31] this leads to transcendental contributions to the anomalous 
dimensions and thus (presumably) to some modification of the underlying Hubbard model 
of [22]. 

An interesting common feature of these observations is the relation of the Hubbard model 
coupling to the Yang-Mills coupling. This raises the possibility that there may be some 
integrable extension of the Hubbard model which contains both elements as part of a larger 
description of iV = 4 super Yang-Mills theory. We will not construct such a model here but we 
will discuss a general approach to constructing a number of supersymmetric Hubbard models. 
Each of these models can be treated perturbatively and thus gives rise to an integrable long- 
ranged spin chain as an effective theory. 

Other supersymmetric generalisations of the Hubbard model have been constructed, see 
e.g. [34,35]. These approaches mainly concern high T c superconductivity models and their 
relation with the t— J model. They essentially use the g/(l|2) or gl(2\2) super algebras, which 
appear as the symmetry algebras of the Hamiltonian of the model. Our approach however is 
different and is based on the QISM framework. It ensures the integrability of the model and 
allows one to obtain local Hubbard-like Hamiltonians for general gl(N\M) superalgebras. 
They can be interpreted in terms of 'electrons' after a Jordan-Wigner transformation. 

The plan of the paper is as follows. In section [21 we define supersymmetric XX models 
whose R-matrices are based on the unitary series gl(N\M). We introduce the corresponding 
Hamiltonians and determine the symmetry of the model. In section [3l we construct the 
associated Hubbard-type model, mimicking the Shastry and Maassarani construction. We 
prove the Yang-Baxter relation for the super Hubbard R-matrix, which allows us to define 
the monodromy and transfer matrices. The symmetry of the super Hubbard model based on 
gl{N\M) is shown to be gl{N -l\M -\)®gl{l\l)®gl{N -l\M -\)®gl(l\l). bisection!! we 
give some examples, writing explicitly the Hamiltonians in the g/(2|2), gl(l\2) and gl(4\4) 
cases. In the first two cases, we also perform a second order perturbation computation 
a la Klein and Seitz [36] and note a relation with the spectrum of the effective two-site 
Hamiltonian with the dilatation operator in the sw(l|2) sector of A = 4 SYM. 

2 Super XX models based on gl(N\M) 

We follow the construction given in [16,21], extending it to the case of superalgebras. In the 
following, we note K = A + M. 

We will use the standard auxiliary space notation, i.e. to any matrix A G End(C K ), we 
associate the matrices A\ = A <g) I and A2 = I <8> A in End(C ) <g> End(C ). More generally, 
when considering equalities in End(C K )® k , we take Aj, j = 1, . . . , k to act trivially in all 
spaces End(C K ), but the j th one. 

To deal with superalgebras, we will also need a Z 2 grading [.] on indices j, such that 
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[j] = will be associated to bosons and [j] = 1 to fermions. Accordingly, the elementary 
matrices (with 1 at position and elsewhere) will have grade [Ey] = [i] + [j], and 
we will use the super-trace: 

K K 

trA = ^(-L)^ if A=J2A ij E ij . (2.1) 

j=l i,j=l 

The grading we use is given by 

b)=\° 1£j£iV ' (2.2) 

11 [1 for iV<j<iV + M. v ! 

2.1 R-matrix 

The R-matrix of the gl(N\M) XX model is defined as: 

R 12 (X) = E 12 P 12 + E 12 sinA + (I <g) I - E 12 ) P 12 cos A (2.3) 
where P\ 2 is the permutation operator, 

K 

Pi2= y 52{-l)®E ij ®E ji (2.4) 
and Si 2 is built from projection operators 

S12 = 7Tl 7T 2 + 7Ti 7T 2 with 7T = ^ £yj , 7T = I - 7T = + . (2.5) 

It is easy to show that E 12 is also a projector, E^ 2 = E 12 . 
Let us introduce the diagonal matrix C: 

C= E jj -E NN -E KK = n-n. (2.6) 

This matrix obeys C 2 = I and is related to the R-matrix through the equalities 

S 13 =^(1-C 1 C7 3 ) and I ® I - E 12 = 1(1 + CiC 3 ) . (2.7) 

One has 

Theorem 2.1 T7ie matrix A2.3\) satisfies the following properties: 

- C-invariance: 

C 1 C 2 R 12 (X) = R 12 (X)C 1 C 2 (2.8) 

- C -parity: 

R 12 (-X)=C 1 R 12 (X)C 2 (2.9) 
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- Symmetry: 

- Unitarity: 

- Regularity : 

- Exchange relation: 



Ri2(\) = P 2 i(A) (2.10) 
Pi 2 (A) R 2 i(-X) = (cos 2 A) I® I (2.11) 
P 12 (0) = P 12 (2.12) 



Rn(X) R2i(fi) = RM #2i(A) (2.13) 

- Yang-Baxter equation (YBE): 

i?i 2 (A 12 ) Pi 3 (A 13 ) P 23 (A 23 ) = P 23 (A 23 ) i?i 3 (Ai 3 ) Pi 2 (A 12 ) where \ tJ = A-A,-. (2.14) 

- Decorated Yang-Baxter equation (dYBE): 

i? 12 (A' 12 ) d P 13 (A 13 ) i? 23 (A 23 ) = i? 23 (A 23 ) i?i 3 (A 13 ) d Pi 2 (A' 12 ) u»& A^- = A, + A,, 

(2.15) 

Proof: C-invariance, C-parity, symmetry, unitarity relation, regularity and exchange relation 
follow from a direct calculation, using the properties 

C 2 = I ; C\ Si 2 = Si 2 C\ = — Si 2 C 2 = — C 2 Si 2 . (2.16) 

The decorated Yang-Baxter equation is a consequence of the Yang-Baxter equation and the 
invariance property. Indeed, the Yang-Baxter equation reads, with the change of variable 
A 2 — > — A 2 , 

Pi2(A' 12 ) i?i 3 (Ai 3 ) -R 23 (— A 23 ) = P 23 (— A 23 ) i?i 3 (Ai 3 ) Pi 2 (A' 12 ) . (2-17) 
Using the antisymmetry property (12.91) . one gets 

i?i 2 (A' 12 ) i?i 3 (A i3 ) C 2 P 23 (A 23 ) C 3 = C 2 P 23 (A' 23 ) C 3 i?i 3 (Ai 3 ) i?i 2 (A' 12 ) . (2.18) 

Multiplying this last equation by C\ C 2 on the left and by C 3 on the right, and using the 
invariance property (12.81) . one obtains (12. 15ft . It remains thus to show the YBE. 
To prove YBE, one evaluates the difference LHS — RHS of f)2.14p . One notes first that 
the terms in S af) , Yj ah P ah and (I <E> I — £ afe )P afe alone satisfy the Yang-Baxter equation. The 
expression is further simplified using the fact that £ a b is a projector and ordering all terms 
with the S's on the left and the P's on the right. One is left, after some algebra and the use 
of standard trigonometric relations, with only two terms: 

Q ; i((Si 2 Si 3 + S 12 S 23 — £ 12 )P 12 P 13 — (S 12 S 23 + S 13 S 23 — S 23 )P 12 P 23 ) (2.19) 
«2(Si 2 — S 23 — S 12 S 13 + S 13 £ 23 )P 13 (2.20) 

where ct\ = sinAi 2 cosAi 3 + cosAi 3 sinA 23 — cos Ai 2 sin Ai 3 cos A 23 and a 2 = sinAi 2 (l — 
cos Ai 3 ) sin A 23 . A direct calculation of these two terms gives identically zero thanks to 
the relation Si 2 = |(1 — CiC 2 ). This ends the proof of YBE. ■ 
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Special case of gl(l\l) 

In the case of gl(\\l) the above construction leads to a trivial R- matrix, because there is no 
index j such that j ^ N, K. However, one can check that modifying the definitions of the 
projectors and C according to 

tt= E n ; 7r = I-n = E 22 ; C = tt-tt (2.21) 

all the properties remain valid. The R- matrix keeps the same form (12. 3ft . with E12 defined 
as in (12.51) . We will use this R-matrix for this particular case. Explicitly, one has 

R(X) = E 2 i <8> E 12 - E 12 ® E 2 i + sin(A) (E n <g> E 22 + E 22 ® E n ) 

+ cos(A) (En ® E u + E 22 <g> £ 22 ) . (2.22) 

2.2 Monodromy and transfer matrices 

From the R-matrix, one constructs the (L sites) monodromy matrix 

£o<i...-L>(A) = -Roi(A) -R 2(A) • • • -Rol(A) (2.23) 

which obeys the relation 

Roo'(X - fj) £ <i..x>(A) £ '<i...l>(aO = £ '<i...l>(aO £o<i...l>(A) i? o'(A ~ A*) • (2.24) 

This relation allows us to construct an {L sites) integrable XX spin chain through the transfer 
matrix 

ti... L {\) = tr £o<i...l>(A) = tr (i? i(A) R 02 (X) • • • Rol(X)) , (2.25) 

where tro is the super-trace in the auxiliary space 0. Indeed, the relation (12.241) implies that 
the transfer matrices for different values of the spectral parameter commute 

[ti... L (X) , ti... L (fi)] = . (2.26) 
Then, the XX-Hamiltonian is defined by 

H = t 1 ... L (0)- 1 t' 1 ... L (0) (2.27) 
where the prime ' denotes the derivative w.r.t. A. Since the R-matrix is regular, H is local: 

L 

H = ^ H jd+1 with H jJ+1 = Pjj+i R' jJ+1 (0) = P jJ+ i Ejj+i (2.28) 
i=i 

where we have used periodic boundary conditions, i.e. identified the site L + 1 with the site 
1. Explicitly, the two-site Hamiltonian reads 

= [ EiN ® Em - EiK ® Eki + ( Em ® Em + Eki ® E ^)} ■ ( 2 - 29 ) 

ijtN,K 
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2.3 Symmetry of super XX models 

Starting from a general K x K matrix A4 generating (a representation of) the superalgebra 
gl(N\M), a direct calculation shows that for 

M = nMn + %M% e gl(N-l\M- 1) © gl(l\l) (2.30) 

we have 

(Mx + M 2 ) R 12 {\) = R 12 (\) (Mi + M 2 ) . (2.31) 

In words, the R-matrix admits a gl(N — 1\M — 1) © gl(l\l) symmetry superalgebra whose 
generators have the form 

E jk ,j,k^N,K for gl(N-l\M-l) 
Enn ; E NK ; E K n ; E K k for <?/(l|l). 

Let us remark en passant that the associated symmetry group is in fact a supergroup, i.e. 
parameters entering the group generators have to be graded according to the grading of the 
superalgebra. This does not affect the 'bosonic' subgroup GL(N — 1) © GL(M — 1) © U(l) © 
£7(1), but the (other) 'fermionic' generators need to have Grassmann valued parameters. 
Note that C-invariance is just a particular case of the above (bosonic) symmetry group. 

As a consequence, the transfer matrix also admits gl(N — 1\M — 1) © gl(l\l) symmetry 
superalgebra, where the generators are given by 



(2.33) 



where M is one of the generators given in (12.321) . The same is true for any Hamiltonian H 
built on the transfer matrix. 

The remaining generators which would allow one to enlarge the symmetry to a gl(N\M) 
superalgebra are given by 

V = ir Mn + ?r Mir generated by E jN ; E jK ; E Nj ; E Kj ,j ^ N,K . (2.34) 

They obey 

VC=-CV so that V x E 12 + E 12 V x = V x ; V 2 S 12 + £ 12 V 2 = V 2 . (2.35) 
This proves that 

{V x + V 2 ) R 12 {\) = R 12 {\) (V, + V 2 ) (2.36) 

where i?i 2 (A) is deduced from Ri 2 (X) by exchanging S 12 and I © I — S 12 . Hence, V is not 
associated to a symmetry of the R-matrix in the usual way. 
Note however that we have the relation 

V x V 2 R 12 (X) = R 12 (X) Vx V 2 . (2.37) 

It induces a gl(N\M) symmetry superalgebra for the XX Hamiltonian, with generators 
Mx ■ M 2 ■ • ■ M.l, where M = Ej%., 1 < j, k < K. Unfortunately, the action of the generators 
on the Hamiltonian eigenvectors is identically zero, except on the pseudo-vacuum. This 
symmetry thus yields no information. 
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2.4 Generalisations 

One can construct a more general R-matrix, defined by 

^i2(A;gi,fe J e 1 ,e 2 ) = £12(51, 52, e u e 2 ) sin A + (e 12 + (I ® I - E 12 ) cos A J P 12 (2.38) 
where 

£12(51, 52, ei, £2) = V] {gi^AfTV ® Ejj + —Ejj ® + q 2 E KK ® + — % ® 

^ 1 ^ 52 J 

+ S { ei ® + — ® ^vtv) + e 2 (q 2 E KK ® % + <g> } (2.39) 



N<j<K 

The parameters (ft, g 2 are complex numbers, while ex, e 2 take values in { — 1,1}. One has 
£ 12 (1,1,1,1) = £ 12 

£12(51, g2,ei,e 2 )Ei 2 (pi,P2, ^1,^2) = Si 2 (gip 1 ,g 2 p 2 ,ei//i,e 2 /i 2 ) . 

Note that only S 12 is a projector. 

It can be checked that the theorem 12.11 is also valid for the R-matrix (I2.38P , except for 
the symmetry (12. 10)) which now reads 

#21 (A; q h q 2 , €l , e 2 ) = i?i 2 (A; — , — , e 1 , e 2 ). (2.40) 

9i 52 

In fact, i?i 2 (A; gi, g 2 , ei, e 2 ) is the (Drinfeld) twist of i?i 2 (A; 1, 1, 1, eie 2 ): 
D 1 R 12 (X; 51, 52, ei, £2) AT 1 = #12 (A; 1, 1, pei, ^e 2 ) 

with D = — Ptvv + — E K k + XI ^ + ^ E H • ( 2 ' 41 ) 

^ ^ 2 j<N N<j<K 

Since D belongs^ to the group SU(N) ® SU(M), the properties proved above (in particular 
the Yang-Baxter equation) remain valid for the matrix i?i 2 (A; gi, g 2 , ei, ei). In the same way, 
it is sufficient to work with the matrix i?i 2 (A; 1, 1, 1, —1) to get the properties of the matrices 
-Ri2(A; gi, 52, ei, — ei), so that there are essentially two different solutions, corresponding to 
the cases e\ = e 2 and e± = — e 2 , hence defining two classes of super XX spin chains. Below, 
we will focus on the R-matrix built on E 12 . 



3 Super-Hubbard models based on gl(N\M) 

We use the R-matrices defined above to build generalisations of the Hubbard model. The 
usual Hubbard model is obtained when we specialise to the case of g/(l|l). We will use the 
results given in [9], generalising them to the case of superalgebras. 

1 

Strictly speaking it is I - ) D which belongs to this group. 
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3.1 R-matrix for super Hubbard models 

One introduces the R-matrix of the super Hubbard model as the coupling of two super XX 
models, according to 

#<i2x34>(Ai, A 2 ) = #13^12) #24^12) + : 7TP\ tanh(/i' 12 ) #i3(A / 12 )Ci#24(A / 12 )C , 2 (3.1) 

sm(A 12 j 

where again A12 = Xi — A 2 and A' 12 = Ai + A 2 . The definition of the parameter h' 12 = 
h(Xi) + h(\2) is given below. It is easy to show that this R-matrix is symmetric 

#<12x34>(Al, A2) = #<34Xl2>(Al, A2) , (3.2) 

regular 

-R<i2><34>(Ai, Ai) = -P<i2><34> = #13 #24 (3-3) 
and obeys the unitarity relation 

#<12><34>(Ai, A 2 )i2<34><i2>(A 2 , A x ) = (cos 4 (Ai 2 ) - f ~t~7T7~7 tanh(/i' 12 )) j I<i 2 > ® I< 3 4> 

V. Vsm(A' 12 ) / J 

where I <13> = I® I. (3.4) 



Property 3.1 When the function h(X) is given by sinh(2/i) = U sin(2A) for some (free) 
parameter U , the R-matrix \3. 1\) obeys YBE: 

#<12x34>(Al, A 2 ) -R<12><56>(Al, A3) -R<34><56> ( A2 , A3) 
— -R<34><56>(A2j A3) -R<12><56>(Al, A3) #<12x34> (Ai,A 2 ). (3.5) 

In that case, the coefficient in jjj3.4\ ) can be rewritten as 
where /i 12 = h(Xi) — h(X 2 ). 

Proof: We use a generalisation to superalgebras of the proof by Shiroishi [37], following 
the proof for algebras presented in [9]. The starting point is the use [38] of the following 
tetrahedral relation [39]: 

l 

#12 #13 #23 = ^2 ^e} #23 #13 #12 > V fl, 6, C = 0, 1 (3.7) 

d,e,f=0 



s 



where R® k = Rjk(Xj — A&), Rj k = Rjk(Xj + A fc ) Cj and the R-matrix is given by (12.31) . The 
non- vanishing entries of the matrix S are given by 

r,0,0,0 _ i 

"0,0,0 ~~ 1 ! 

q1,1,0 _ -| q0,1,1 _ i r-i.O- 1 _ i 

"1,1,0 ~~ 1 ' 0,1,1 ~~ 1 i "1,0,1 ~~ 1 

-,1,0,0 _ t//\ \ \ \ chOfi _ u//\ \ \ \ cl,0,0 



^o,'o'i = ^(Ai,A 2) — A 3 ) ; S y = W(Xi, A 2 , — A3) ; S^'i = U(\\, A 2 , — A3) 

^o,'o,'i = ^(Ai, — A 2 , A 3 ) ; ^'q'q = V^(Ai, — A 2 , A 3 ) ; Si^'i = V(Ai, — A 2 , A 3 ) 

^o.'i'o = U(— Ai, A 2 , A 3 ) ; S'ioo = V(— Ai, A 2 , A 3 ) ; 5 r °i : 1 = W(— Ai, A 2 , A3) 

^o.'o.'i = ^(Aij A 2 , A 3 ) ; ^l.'o = ^(Ai, A 2 , A 3 ) ; ^1,0,0 = ^(Ai, A 2 , A3) 



(3. 



with 



U{\i, A 2 , A 3 ) = - C °fj Sin if 23 | ; V{\ 1 , A 2 , A 3 ) - Sin(A ' 12) Sin(A23: 



sin(Ai 3 ) cos(A 23 ) ' sin(Ai 2 ) sin(A 



23 , 



W{K A2 ' As) = rin(Aij cos(Alj ; X]k = ^ ~ h ! X * k = Xj + X ^ 3, k = 1,2,3 . (3.9) 
One needs also the relations: 

pi p 1 p 1 - sin ( A i3) cos (A 23 ) o ! cos(A / 12 ) cos(A 23 ) x 

-"■23 - ft 13 /t 12 — „;„/\ N "\ 23 -^13 12 "i „_„/ \ \„^„/\ \ -"-23 -"-13 /t 12 

COS(Ai2j COS(A23j 

^23^13^12! (3.10) 

i P 1 P 1 R° _i_ cos (Al2) CQS(A 23 ) x j 
^ ^23 «13 ^12 + cos(A / a) ^23 ^13 ^12 

^23^13^12- (3-11) 



sin(Ai 3 ) 


cos(A 23 ) 


cos(A' 12 ) 


sin(A' 13 ) 


' cos(Ai 2 ) 


sin(Ai 3 ) 


sin(Ai 3 ) 


cos(A 23 ) 


sin(A' 13 ) 


cos(A' 23 ) 


cos(A i2 ) 


sin(Ai 3 ) 



cos(A' 12 ) sin(Aj 



13; 



It is easy to prove, using for instance a symbolic computer program [40], that all these 
relations hold for the R-matrix (12.31) . provided Cj = 1 and P\i is a (super) permutation 
operator. 

Then, the end of the proof is similar to the algebra case: a direct (but lengthy) calculation 
shows that the matrix 

-R<i2><34>(Ai, A 2 ) = -Ri3(A i2 ) i? 2 4(Ai 2 ) + a(Ai, A 2 ) i2 13 (A' 12 ) i? 24 (A' 12 ) C\ C 2 (3.12) 

obeys YBE provided the matrix i? 12 (A) obeys theorem 12.11 relations (13. 7p and fl3.10H3.lT]) . 
and a(Ai, A 2 ) is given by 

cos(Ax - A 2 ) sinh(/ii - h 2 ) , . 

«(Ai,A 2 ) = tt j—r -— — , (3.13) 

v ' cos(Ai + A 2 ) cosh(/i! + h 2 ) 

where hj = h(Xj), j — 1,2 is defined by sinh(2/i) = U sin(2A) ; see [9] for more details. ■ 
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3.2 Monodromy matrices, transfer matrices and Hamiltonians 

We remind the reader of the usual proof of integrability for models based on transfer matrices. 
Let TZ a b(Xi, A 2 ) be an R-matrix obeying YBE, and being regular (7Z a b(X, A) = P a b)- a and 
b denote the 'coupled' spaces (a, b =< 12 >, < 34 > in the above cases). From YBE, one 
deduces that the monodromy matrix 

£a<bi...b L > (Ai,A 2 ) =K abl (\ 1 ,\ 2 )...K abli (\ 1 ,\ 2 ) (3.14) 

obeys 

7^aa'(Ai, A 2 )£ a (Ai, A 3 )£ a /(A 2 , A 3 ) = £ a /(A 2 , A 3 )£ a (Ai, A 3 )7?. aa /(Ai, A 2 ) (3.15) 

where the dependence in the quantum spaces b\, . . . , b^ has been omitted in C This relation 
proves that one can define a transfer matrix 

t(A l5 A 3 ) =tr a £ Q (A 1 ,A 3 ) (3.16) 

which obeys 

[t(Ai,A 3 ) J(A 2 ,A 3 )] =0. (3.17) 
From the transfer matrix, one then deduces that all the Hamiltonians 



JT( A *)=f(0 ) /*)- 1 ^*(A,/*) 



(3.18) 



A=0 



define, for any A, an integrable model, since we have 

[H(n),t(\,(jL)] = 0, VA. (3.19) 

However, demanding further that the Hamiltonian be local, one is led (using the regularity 
property) to specify n = 0. One then gets 

[H, t(X)} = , VA , for H = H(0) = t(O)" 1 £'(0) and t(X) = t(X, 0) . 

The transfer matrix t(X) is constructed from the 'reduced' monodromy matrix 

L a<bu . ML> (X) = TZ abl (X,0)...TZ abL (X,0). (3.20) 

This 'reduced' monodromy matrix is just the one used to define the Hubbard model; one 
can compute 

J R<i2><34>(A,0) = J R 13 (A) J R 24 (A)(l®I + tanh(/ i )C 1 C 2 ) . (3.21) 

Hence, it is the locality requirement that imposes the form of the monodromy matrix used 
for the Hubbard model. More general (a priori non local) Hamiltonians can be defined using 
the form fl3TT8l . 

The matricial form of the Hubbard-type Hamiltonian reads 

L L 
H = ^2 H<2j-l,2jX2j+l,2j+2> = ^2 ^<2i-l,2j><2i+l,2j+2>(0) (3.22) 

i=i i=i 
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with R ab = P ab R ab , R' ab (0) = ^Rab(\ 0) 



and 

A=0 



H<2j-l,2j><2j+l > 2j+2> = E 2 j-l,2j+lP2j-l,2j+l + ^2 j,2j+2 P2j,2j+2 + — (C 2 j-lC 2 j + C 2 j+lC2j+2^ 

(3.23) 

where we have used periodic boundary conditions, i.e. identified sites < 2L + 1,2L + 2 > 
with sites < 1,2 >. 

Remark 3.1 Let us remark that, due to the coupling of the two XX models, the total 
number of sites is 2L, but the number of 'coupled' sites (which are the real physical ones) is 
L. We will thus refer to the Hubbard Hamiltonian f)3.22p as an L-site Hamiltonian. This is 
consistent with the notation used after Jordan- Wigner transformation (see below). 

3.2.1 Gauged version of the super Hubbard model 

In the literature [12, 17, 18], a gauged version of the Hubbard R-matrix is used. It is defined 
by 

Dg ( \ \ \ _ AhidC'2 e \h 2 C 3 C A p (\ \ \ -5/11C1C2 e -hh 2 C 3 C4, 

- K <i2><34>l A i> A 2) — e 2 e 2 -n:<i2><34>lAi, A 2 J e 2 e 2 

where hj = h(\j) , j = 1,2 (3.24) 

By construction, -R<i 2 ><34>(^i! ^2) also obeys YBE, and is unitary, symmetric and regular. 
Following the same steps as before, we introduce the 'reduced' R-matrix 

^ S <i2><34>(A,0) = —1^ I 12 (h) R 13 (\) R 2i (X) I 12 (h) , (3.25) 

where 

h h 
I 12 (h) = cosh(-) I ® I + sinh(-) d C 2 . (3.26) 

It leads to the same Hamiltonian fl3.22l) - fl3.23l) . This gauged version was originally introduced 
to recover the exact form of Shastry's R-matrix. 

3.3 Symmetries 

We generalise to superalgebras the results obtained for su(N) Hubbard models (see for 
instance [9,17]). For completeness, we compare them with the well-known symmetry of the 
usual Hubbard model [7,37]. 

Proposition 3.2 The transfer matrix of the Hubbard model admits a gl(N — 1\M — 1) © 
gl(l\l)@gl(N-l\M-l)@gl(l\l) symmetry algebra, each of the gl(N - 1\M - 1) ©flfZ(l|l) 
corresponding to the symmetry of one XX model. 

As a consequence this symmetry is also valid for the Hubbard Hamiltonian. 

Proof: To prove this symmetry, it is sufficient to remark that 

MC = CM (3.27) 
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where M is given in (12.301) . Thus, one gets 



(A, 0) , Mi + M 3 ] = = [ J R<i2><34> (A, 0) , M 2 + M 4 ] (3.28) 

where i?<i2><34>(A, 0) is the R-matrix of the Hubbard model. 

As far as Hamiltonians are concerned, the generators of the symmetry have the form 

L L 

M evn = ^M 2j and M odd = ^M 2j _i (3.29) 

3=1 3=1 

They generate a gl(N - 1\M - 1) © gl(l\l) © gl(N - 1\M - 1) © gl(l\l) superalgebra. ■ 

It is well-known that the Hubbard model possesses a gl(2) © gl(2), and thus it is natural 
to look for a gl(N\M) © gl(N\M) symmetry algebra for the generalised Hubbard models. 
Unfortunately, it seems not to be present. To discuss this point, we now review how the 
so(4) symmetry algebra is obtained in the framework of the Hubbard model and point out 
some properties which are valid only in this case. 

3.3.1 Enhancement of the symmetry for Hubbard model 

As has been shown (originally in [7], see also [9]), the full symmetry of the periodic Hubbard 
model (for finite L) can be obtained through a change of the Z 2 -grading. In the present 
context, it amounts to consider the g/(l|l) superalgebra[f|: 

Proposition 3.3 In the gl(l\l) case, the Hubbard R-matrix obeys 

(Vif — V^) -R<i2><34>(Ai, A 2 ) = — -R<i2><34>(Ai, A 2 ) (V 12 ~ ^34) (3.30) 
(W± + W£) #<i2><34>(Ai, A 2 ) = i?<i 2 x34>(Ai, A 2 ) (W^ + W£) (3.31) 

where V ± = a±® a± and W ± = a±® <j t . 

These relations are not valid any more for a general gl(N\M) superalgebra for generators 
V = E Nj © E N j, E K j © E K j, E jN © E jN or E jK © E jK and W = E Nj © E jN , E Kj © E jK , 
E jN © E Nj or Ej K © E Kj . 

Proof: Direct calculation. In particular, we checked that this relation does not hold for 
gl(l\2). m 

Relations (13.30p — (13.31)) are then enough to deduce the following corollary, proved in [9]: 
Corollary 3.4 For gl(l\l), the Hamiltonian 

L 

H = / ] -^<2.7-l,2,7><23+l,2i+2>(Q) (3.32) 
3=1 

possesses a gl(2) © gl(l) © gl(l) symmetry algebra when L is odd; this symmetry extends to 
a gl{2) © gl{2) algebra when L is even. 

3 It corresponds to the R matrix in the notation of [9], 
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The generators of this symmetry have the form 

L L 

= E(-^ and 5f° = E w S-w ^ 3 - 33 ) 

3=1 3=1 
2L L 

5? = E C 'i anrf ^^E^i-^)- (3.34) 
3.3.2 Comparison with the l gl(2) Hubbard model' 

We give here the counterpart of the sect ion 13 . 3 . 1 1 when dealing with the gl(2) Hubbard model, 
constructed using the transfer matrix approach. 

Proposition 3.5 In the gl(2) case, the Hubbard R-matrix obeys 

(^12 C3C4 — V34) -R<12><34>(Al, A2) = — P<12X34> (Ai,A 2 )(^-C7iC 2 ^) (3.35) 
<12><34>(Al, A 2 ) — -R<12><34> (Ai, A 2 ) {W? 2 + C 1 C 2 W&) (3.36) 

where V ± = a±® a± and = a±® a T . 

These relations are not valid any more for a general gl(N) algebra for generators V = 
E Nj ® E N j or E jN (g) Ej N and W = E Nj ® E jN or E jN (g) E Nj . 

Proof: A direct calculation shows that the relations 

7Tl ^2^13^4 = ^3^4^24 ^4 ) ^^P^V^ = ^1 ^2 ^24 V 2 % ] 7Tl7f 2 P 23 V 2 \ = TT 3 TTiP U V^ 

n^P^W^ = H^PiaWu ; 7ri7r 2 Fi3Wi4 = ^ir^P^W^ ; 7ri7r 2 Pi 3 W^ = ir^P^W^ 

hold for gl{2), but not for the other (super) algebras. Using these relations, it is then easy to 
deduce the relations H3.35f) — ( T3.36D . We also checked by direct calculation that the relations 
- f[33ol) do not hold for gl(3). m 



Corollary 3.6 For gl{2), the Hamiltonians 

L-l 

(Ai,A 2 ) (3.37) 

3=1 

have a gl(2) © gl{2) symmetry algebra. 

R implies the same symmetry for the non-periodic Hubbard Hamiltonian 

L-l 

Hn.p. — E -^'<2i-l,2jX2j+l,2j+2>(0) ■ (3.38) 
3=1 

Proof: Multiplying from the left by P i3 P 2 4, the relations fl3.35p - fl3.36p can be recast as 

(Yl2 ~ ^34 CiC 2 ) P<12><34>(Al, A 2 ) = P<12><34> (X l ,\2)(V±-C 1 C 2 V±), (3.39) 
(W± + W£ CiC 2 ) P<i 2 ><34>(Ai, A 2 ) = R <12><34> (Ai, A 2 ) (Wf 2 + CtC 2 W£) (3.40) 
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It shows that the generators (again with V = a±® a± and W = a±® a T ) 

L 

= V£ -C 1 C 2 V± + C 1 C 2 C 3 C 4 V±-C 1 C 2 C 3 C 4 C 5 C 6 V T i + ... (3.41) 

= W± + CiC 2 w£ + CiCsCs^ W^ 6 + C&CsCtGiCe W& + ... (3.42) 

commute with the above Hamiltonian (with no restriction on the parity of L). It is trivial 
to check that they form a gl(2) © gl{2) algebra, with Cartan generators 

L 

= X)(qw-i ±C ^)- (3.43) 



The Hamiltonians H(\i,X 2 ) and i? n .p. are not periodic, since they do not contain the 
term -R<2l-i,2L><i,2>(Ai, A 2 ), which breaks the symmetry. Hence, H np _ does not correspond 
to the usual Hubbard model. However, in the thermodynamical limit L — > 00, the missing 
periodic term is sent to infinity, and one recovers the symmetry of the usual Hubbard model. 



3.3.3 Jordan— Wigner transformation and periodicity 

Anticipating the reminder of section 13.51 on Jordan-Wigner transformation [41], one is 
tempted to associate the gl(2) construction of [12, 14] to the Hubbard model, but it is 
well-known that, for algebras, the Jordan-Wigner transformation does not preserve the pe- 
riodic boundary condition [13] (see also [9]). Indeed, through this transformation, one gets 
for instance 

c) c j+1 - E® Eg +1) , j = 1, 2, . . . (3.44) 

where the superscript indicates the site to which the matrices belong, and the arrow denotes 
the Jordan-Wigner transformation. From periodicity, one should thus get 

clc^E^E^ (3.45) 

However, performing the Jordan-Wigner transformation, one gets 

clc 1 ^E$E$(C 1 ---C L _ 1 ) (3.46) 

Hence, in the gl(2) case, the Hubbard Hamiltonian we obtain is non-periodic in terms of c 
and (i.e. after Jordan-Wigner transformation). 

When dealing with superalgebras, the Jordan-Wigner transformation is modified [9] (see 
a reminder in section l3~5|) . and now respects the periodic boundary condition. Due to this, 
we obtain the usual (periodic) Hubbard Hamiltonian in the case of gl(l\l). 
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In other words, for algebras, the Jordan-Wigner transformation needs to modify the 
bosonic/fermionic character of some operators: this is done using (non-local) products of 
Cj = (1 — 2rij) generators which break the periodicity. For superalgebras, no change of char- 
acter is needed; the transformation is a local isomorphism, so that periodicity is preserved. 
In this respect, the superalgebra case looks more natural than the algebraic one. 

These considerations are consistent with the results of sections 13.3. II and 13.3.21 about the 
symmetry of non-periodic gl{2) and periodic gl(l\l) Hubbard models. 



3.4 Change of notation 

The above presentation of the Hubbard model is based on the transfer matrix formalism, 
the Hubbard model itself being obtained by coupling two independent XX models, hence 
the notation used for the Hubbard Hamiltonian (13.221) . In the following, we are dealing 
with explicit expressions of this Hamiltonian in specific cases and we would like to make 
contact with the notation commonly used in particular in the condensed matter community. 
Therefore, we will perform a change of notation in the rest of the paper in order to stick to 
more familiar expressions. 

The construction of the Hubbard Hamiltonian, see eqs. fl3.22p - fl3.23l) . shows that one 
considers a 2L site lattice on which live two independent XX models, the first one living on 
the odd sites, the second one on the even sites. We introduce a map on the site labels in 
such a way that the 2L site lattice of the coupled XX models is interpreted as a L site lattice 
for the Hubbard model: 

<2j-l,2j> - (J = 1,...,L) (3.47) 

the operators living on the odd (even) sublattice being labelled by j (J,). With this notation, 
the Hubbard Hamiltonian fl3~2"2]) - fl3~2^|) reads 

i=i 

(3.48) 

where we used the periodicity conditions. 



3.5 Jordan— Wigner transformation 

Let us consider p sets of fermionic oscillators cf \ cf (i = 1, . . . , L and q = 1, . . . ,p) that 
satisfy the usual anticommutation relations 

{cf\ cf = 5 l3 5 qql { C W cf } = c ft } = o (3 . 49) 

One defines the following matrix (where the usual number operator) 

X i q) = ( 1 c («)t ^9) J ' ( 3 - 50 ) 
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The entries xj?^a °f this matrix have a natural grading given by [a] + [/?] where [1] = 1 and 
[2] = 0. 

In the gl(2 p 1 \2 P x ) case, one defines at each site i the generators 

p a— 1 

-V-:,,.,.,../..,.;. = (-IT X^ ia , . . . X^ pa , where s = j> fl ] ( £ (H + K])) . (3.51) 

a=2 6=1 

It is easy to verify the following properties: 

(^i;Qi...Q p ,o' 1 ...Op) -Xi;a' 1 ...a' p ,ai...a p (3.52) 

^i;ai...c*p,a' 1 ...a£ ^i;/3i.../3 p ,^...^ = • • • <^/3 p X; Ql ...a p ,/^.../^ (3.53) 

ap,ai... a p = 1 (3.54) 

ai,...,a p 

■^t;ai...a p ,a / 1 ...c4 -Xj;/3i.../3 p ,#.../% = (-l) 5 ^';/?!...^,^...^ ^;ai...a p ,ai...a p (« 7^ j) (3.55) 



where p = ([aj + k])) ([ft] + [ft 

o=l 6=1 



The first three properties are local (on site) while the last one relates different sites. They 
state that the X ;a ' 1 ,..., a j ) ,o 1 ,...,o J , form an algebra isomorphic to the tensor product of L copies 
of gl(2 p ~ l \2 p ~ 1 ) . The mapping X — > E is known as a Jordan-Wigner transformation. We 
observe that it can be uniquely defined once an entire line is given; indeed, hermitian 
conjugation fixes the corresponding column so the full matrix can be reconstructed in the 
following steps: 

(i) 

1. the element E\{ is associated to one of the 2 P diagonal generators -X» ;ai) ... ja aij ... )a ; 

2. the remaining 2 P ~ 1 — 1 bosonic generators are freely associated to the bosonic ones 
E?l a = 2,...,2T\ 

3. the 2 P ~ 1 fermionic generators are freely associated to the E^, a = 2 P ~ 1 + 1, . . . , 2 P . 

All specific realisations are isomorphic because they can be obtained one from the other by 
exchanging lines and columns of the matrices. An example of such a mapping is given in 
( 14. 2p . The gl(N\M) cases that are not of the form gl(2 p ~ 1 \2 p ~ 1 ) are "incomplete" and can 
be obtained by embedding in the smallest algebra gl(2 p ~ 1 \2 p ~ 1 ) such that N < 2 P ~ 1 and 
M < 2 P ~ 1 . Then, by removing an appropriate choice of lines and columns, one projects the 
matrix Xi ]0 ,' it ___ t0 ,^ ait ___ j0tp to a gl(N\M) subalgebra. This can be done in many inequivalent 
ways and has been done in section fl~2l with the projector ( 14.24ft . In this sense, any gl(N\M) 
Hamiltonian describes a sector contained in the larger gl(2 p ~ 1 \2 p ~ 1 ) Hamiltonian's space of 
states. 

4 Examples 

4.1 gl(2\2) Hamiltonians 

In the gl{2\2) case, the generators X^ a j3 j0l ipi at each site i are given by 

X a ,ai — (— 1 ^H+KDM Y , X' (4 11 
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They are mapped on the E pq matrices with the following assignment of indices (a, (3, a', (3' = 
1,2 and p,q = 1 . . .4): 

11 -> 1, 12 -> 3, 21 -> 4, 22 -> 2 (4.2) 

which respects the grading in the sense that if (a/3, a'j3') — > (p, q), the grades of X a p j0l tpi and 
of E pq coincide. 

Then the gl(2\2) XX Hamiltonian (12.281) reads as (the subscripts correspond to the site 
indices) : 

HS 2) = E (c]c i+1 + cl +lCl ) (d^ i+1 + JlA + 1 - nj - n; +1 ) (4.3) 
i=i 

L 

Ef_ 't t ' 't t ' _i_ t 't ' , 't t ' 

1 C i C i C i+l C i+l C i+l C i+l C i C i + C i C i+l C i C «+l + C j C i+l C « C i+l 

i=l 

+ (1 - nj - n' i+1 ) (cjc t+1 + c[ +lCi ) } . (4.4) 

This Hamiltonian exhibits interesting features. First of all, the number of pairs (i.e. doubly 
occupied sites with one unprimed and one primed particle) is conserved by the Hamiltonian, 
so that one can restrict the study to sectors with a given number of pairs. The first two 
terms of (14.41) correspond to a BCS-like conductivity in the physical space (pair hopping), 
while the last term corresponds to ordinary conductivity (hopping for unprimed particles 
with interaction with a background of primed particles). The middle term corresponds to 
an exchange between the two types of particles. 

As explained in section El the Hubbard-type Hamiltonian (13.221) is obtained by coupling 
two copies of XX Hamiltonians, with fermionic oscillators c^ ai and c a ^, a =f , |. Hence, one 
gets for the g/(2|2) Hubbard Hamiltonian: 

L 

HfluP = E { E ( C i,i Ca ' i + 1 + C i,i+l C °,i) { C a,i C 'a,i+l + C 'a,i+l C 'a,i + 1 _ n 'a,i ~ n 'o,i+l) 
i=l <r=U 

+ Z7(l-2n T , i )(l-2r H , i ) )• (4.5) 



The space of states at each site i is spanned by the vacuum |0)j, the up states | f)j, I Y)u 
1 the down states | | J/)*, | and by tensoring the up states with the down states, 
where \a) t = cj. )< |0) j , \a')i = (Q ti \0)i and \aa') t = cJ^JO^ with a =f, |. 
This Hamiltonian can be compared with the #/(4) Hubbard Hamiltonian which is given by 



H 9m _ V / V (J r • , r n r' 4- r 1 " 

11 Hub ~ \ Z-c \ C o-,i c v,r+l C a,i C a,i+l ~T~ C o-,i 

i=l (T=t,J. 

+ n CT,i n CT,i+l( C l,i C o-,i+l + c l,i+l c <r,i) + n cr,j^o-,i+l(cJ j jC^ i+1 + cJ ji+1 C^j)) 

+ t/(l-2n Ti y Ti J(l-2n i)i n' w ) }, (4.6) 
which is free of exchange terms. 
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It is of interest to make a perturbative calculation of the g/(2|2) Hubbard Hamiltonian 
( 14. 5p a la Klein and Seitz [36]. To this aim, one introduces the notation 



<T=T.J. 



(4.7) 



where M' a>ij = c^c^ + cl/ ai + 1 



n' 



n' a y The Hamiltonian takes then the form 



HfT = + X i+lii ) + UJ2(1- 2n u )(l - 2n L 



i=l 



i=l 



At large U, the potential term is the dominant one, while the X term can be treated as a 
perturbation. From the form of the potential term, one is led to define a projector n on 
singly occupied states with unprimed particles (i.e. | f ) or | j)), without any limitations on 
the primed particles: 



n i,i) 



(4.9) 



i=l 



Then, one can easily check that X\- = Xji and that Ho fulfills the following conditions: 

n x„n = o, (1-^)^x^0 = 0, u x iA+l x i+w u = o . (4.io) 

Note that the 'dressing' factors M' ai j play no role in the derivation of these relations since Ho 
and M' ai j commute. It follows that the effective Hamiltonian at second order of perturbation 
(the first order is vanishing) is given by 



1 L 

i=i 



(4.11) 



After some simple algebra, one finally gets 

L 



H 



(2) 
eft' 



£5> 

i=l 



2S*S? +1 ) ($i + $i ^i+l) ^,11+1^1,11+1 



(4.12) 



C L C t,* and S i = U n ^i- n i,i)- 



where one has defined the s/(2) generators 

It is worthwhile to emphasise that the symmetry algebra of the Hubbard Hamiltonian 
(14. 5 p is g/(l|l)©g/(l|l)©(7Z(l|l)©(7/(l|l), the symmetry algebra of the effective Hamiltonian 
at second order of perturbation (I4.12p however is enhanced to g/(2|2) © g/(2|2). Indeed, the 
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following generators 

L L L L 

Y n ^ c 'h c 'n ' P° = Y C 'i C T,i ' P° = I X n ^K,* - > ^ = X! n ->' 



Pa 



i=l i=l i=l i=l 

L L L 



Va 



Y U °,i C 'li C 'li ' ^ = S ™<M 4* C T,i ' ^ = \ Y n crA n \,i + n U - X ) 
i=l j=l i=l 

L L 

Y U '-T,i C r,i ' <Kt = Y ^ U '-T,i C 'r,i ( r = t, I) 

i=\ i=l 
L L 

Y 71 ^ 1 ~ n '-T,i) C T,i » Xar = XX,^ 1 ~ ™-T,i) C T,i (T=T,J.) ( 4 -!3) 

i=l i=l 

which generate two commuting copies of s/(2|2) (one for a =| and one for cr =|), commute 
with the effective Hamiltonian (14.121) . Due to the presence of the projector Ho in H^ s , the 
number operator n CT| j can be expressed as h±S£ (+ for j and — for |). The generators p^, 
and 77^,77* generate the four commuting sl(2) algebras, and the remaining non-vanishing 
commutation relations are given by 



X.OT 



\pU 


£] = 




[^Xar] = 




(4.14) 


\pU 


£] = 




[Pa.Xarl = 




(4.15) 


k T ,< 


&] = 




K T ,^] = 




(4.16) 




&] = 




feXoV] = 




(4.17) 




>X^} 


= -{^. = 


{4>%,x%} 




(4.18) 


{fa 






\X(JT ' X(JT J 


= \N a - V z a + rpl 


(4.19) 



where a, r =f , j. and e = ±. 

One has to add two more generators (the U(l) factors such that one gets two gl(2\2) super- 
algebras), which are both represented in the present realisation by matrices proportional to 
the unit matrix. 



4.1.1 Spectrum of the Hamiltonian and comparison with J\f = 4-SYM 

It is interesting to check if (14.121) has some relation with the dilatation operator of some 
sector in M = 4-SYM. The proper candidate is su(l|2) sector whose two-site Hamiltonian is 
Hsym = 1 — P12 [42]. On each site, after the projection Ilo, the two-site Hamiltonian (14.121) 
acts on eight states 



IT), 


im, 


im, 


irm 


(4.20) 


II), 


im. 


nn, 


IITT) 


(4.21) 



so that it is a 64 x 64 matrix. It has 32 vanishing lines and columns and the remaining part 
is built of two-by-two blocks of the form 

B - = (-i "/) 01 B * = (I 1) ' (4 22) 
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This means that the eigenvalue zero is represented 48 times and the eigenvalue 2 appears 16 
times. 

The H$ym Hamiltonian is a 9 x 9 matrix with two empty lines and columns, three 
blocks B- and a one-dimensional diagonal entry with value 2. Therefore, our Hamiltonian 
(I4.12p contains the correct su{\\2) spectrum. The interpretation of states is not obvious 
because the one-dimensional block with value 2 is absent and we can obtain it only after 
a diagonalisation of one of the blocks (14.221) namely by mixing the states on two sites. 
Moreover, the enhancement of symmetry seems to be strictly a feature of this second order 
Hamiltonian and is most probably lost at higher orders. 



4.2 gl(l\2) Hamiltonians 

Following formula (12.281) . the g/(l|2) XX-Hamiltonian can be obtained from the gl{2\2) one 
by suppressing the index 1 for example (or equivalently the index 3). One gets therefore 

L 

H xx 2) = { c i c l+i c i c 'i+i + 4 c i+i c 'iCi+i ~ (cjci+i + 4 +1 Ci)« +1 | (4.23) 



where at each site i the space of states is spanned by c]|0), c} |0) and cjcf |0). In other words, 
the space of states of the gl(l\2) XX model can be obtained from the space of states of the 
gl(2\2) XX model by acting with the projector 

L 

= l[(rii +<- m<) (4-24) 

i=l 

It can be easily verified that one has 

tt91(1\2) _ n rrfl«(2|2) n / , r,^ 

H XX - i V(l|2) H XX n s«(l|2) I 4 - 25 ) 

As illustrated in the previous example, the gl(l\2) Hubbard Hamiltonian is constructed 
by coupling two copies of the XX Hamiltonian, with fermionic oscillators c* ai and c a<i , a =f, J,. 
It reads therefore 



H-Hub ~ { ( C i,i C il,i+l Co ".* C cr,i+l C l;i C l,i+l C a,i C o;i+l ( C l,i C cr,i+l + C l,i+l C cr,i) n a ,i n a 

+ U {n' hi - n tl - n r>i n' hi ) (n' u - n Ui - n u n' Ui ) } (4.26) 

the space of states at each site i being spanned by tensoring the up states | | t')i) I TT')i with 
the down states | | \ 11%, where \a) t = cj. ji |0) i , \<r')i = cJjO)^ and \acr , ) i = 4,i c Jil )i 
with cr =|, |. It follows that the space of states of the gl(l\2) Hubbard model is obtained 
from the space of states of the gl(2\2) Hubbard model by acting with the projector 

L 

n 9 i(i|2) = \\i n hi + n \,i - n ^i n \,i) ( n i,i + n \,i - n i,i n 'i,i) ( 4 - 27 ) 

i=l 
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Again, one has 

<f = rW, 2) <f ) n s , (1|2) (4.28) 

which is a direct consequence of (I4.25P and the trivial embedding of the gl(l\2) and gl(2\2) 
C matrices entering the definition of the potential term. 

Introducing the notation = £)o-=ti c i,i c <r,j N' a ,iji see ec L- (14.71) . where now M' a ^ = %^ a ,j ~ 
n' a i n' a j, one has 

L 

H 9 hT = + X»i,i) + ^("r,i - *t.i - n u n' u )(n' Ui - n u - n^n'J (4.29) 

i=l 

At site i, among the nine possible states, four of them have an interaction energy —U and 
the other five have an interaction energy +U. These four states are characterised by the 
constraint n^j + n^j = 1, hence the projector n in the effective Hamiltonian is again given 
by (14.91) and the relations (I4.10p are still satisfied. Therefore the effective Hamiltonian at 
second order of perturbation reads 



n n 9Z(1 | 2) (4.30) 

where the 'dressing' factor M' a ^ is obtained from N' ai i by action of the projector n ff z(i| 2 ): 
M' .. — c t d ■ + c /f d ■ - n' n' ■ 

Unfortunately, the symmetry of this Hamiltonian is not gl(l\2) © gl(l\2) as one might 
have hoped from the gl (2 1 2) case, but only gl(l\\)®U(l)(Bgl{l\l)®U{\). The corresponding 
bosonic generators are given by p*, r]*, N a (with a =|, j) and the fermionic ones by (f)++ and 

The spectrum of this Hamiltonian is completely contained in the g/(2|2) case (14.121) 
already described in the previous section. 

4.3 g/(4|4) Hamiltonians 

In the #Z(4|4) case, the generators X^a^^pi^ at each site i are given by 

The Xa^^ipry generators are mapped on the E pq matrices with the following assignment of 
indices (a, 0, 7, a', 0', 7' = 1, 2 and p,q = 1 ... 8): 

111 1, 112 -> 6, 121 -> 7, 122 -> 4, 211 -> 5, 212 -> 2, 221 -> 3, 222 -> 8 (4.32) 

Again, this assignment respects the grading in the sense that if (a(3j,a'P'j') —> (p,q), 
X a py ta ipiy and -Epg have the same grade. 
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i=l 
L 



Then the gl(4\4) XX Hamiltonian (12.281) reads as (the subscripts correspond to the site 
indices) : 

L 

H gl(4\4) = [c\d +l + c\ +l d + 1 - Ui - U l+l ) (c'?c' i+1 C?c'l +1 + C^c"^ 

i=l 

~ « + i(c!^Ui + - <^lM 4 + i + clA)) (4-33) 

As in the g/(2|2) case, one gets for the g/(4|4) Hubbard Hamiltonian: 

L 

^Hub ^ = { i C lr,i C v,i+l + C lr,i+l C o;i + 1 — n a,i ~ ( K C l,i C a,i+l C J,i C a,i+l 

_i_ 't I "t J' 1 l I "t " I H\ II \ 

' C a,i+l C a,i C a,i+l C a,i n a,i n <j,i+l\ C a,i C <J,i+l ' C a,i+l C a,i) 

~ < 4 <m(<*<m + + U(l - 2n\^{l - 2^X1) } ( 4 - 34 ) 

One observes that this Hamiltonian exhibits a 'Russian doll' structure. Indeed, there are 
four sectors in the space of states where the gZ(4|4) Hamiltonian reduces to the g/(2|2) one. 
These sectors are defined respectively by = n'^ = 1, = n"^ = 1, n"^ = n'^ = 1, 
n '{i = "^'/j = 1 for 1 < z < L. The obtained Hamiltonian can be further reduced to gl(l\l) 
Hamiltonian by imposing on each site n a ^ = or ?v,i = 1- 

5 Conclusion and perspectives 

We have constructed super- Hubbard models based on the superalgebras gl(N\M), with a 
special focus on models that may apply to SYM theories. We have seen that in the case 
of superalgebras, the Jordan- Wigner transformation is a local isomorphism. Therefore, the 
interpretation of the models in terms of 'electrons' is more natural. 

The symmetry superalgebra and the Hamiltonian have been given, and we performed a 
perturbative calculation a la Klein and Seitz [36] for the Hamiltonians based on the super- 
algebras gl{l\2) and gl(2\2). 

The next step in the study of our models is the determination of the spectrum and 
the Bethe equations, as they were constructed for Hubbard or generalisation, using the 
algebraic Bethe ansatz [19-21,43]. This is an heavy calculation which we postpone for 
further publication, but from the analytical Bethe ansatz approach, one can guess their 
form. In particular, as for spin chain models, one expects as many presentations of the Bethe 
equations as there are inequivalent Dynkin diagrams. All these presentations should lead to 
the same spectrum. For more informations, we refer to [44,45] where similar calculations 
were performed in the case of XXX super spin chains. 
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